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Abstract. We improve the Winitzki's Approximation of the error function 

I 2 A-|-0.147x 2 

erf(x) = V 1 - e~ x 1+0.147* 2 which has error \e{x)\ < 1.25 • 10~ 4 Vx > till 
reaching 4 decimals of precision with |e(x)| < 2.27 • 10 -5 ; also reducing slightly 
the relative error. Old formula and ours are both explicitly invertible, essentially 
solving a biquadratic equation, after obvious substitutions. Then we derive 
approximations to 4 decimals of normal cumulative distribution function 
of erfc(x) and of the Q function (or cPhi). 
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In this note we improve the Winitzki's Approximation of the error function till 
reaching 4 decimals of precision both for erf and reducing about 5.5 times 
the respective absolute errors (and reducing the relative errors too). 
Lemma (Winitzki's Approximation of erf). (See [1]) 



erf(x) = V 1 - e x 1+0.147** < 1.25-10~ 4 |e r (i)| < 1.28- 1CT 4 Vx > 

Theorem (Improving of Winitzki's Approximation of erf). 



-1.2735457* 2 -0.1487936x 4 
BVf(x) = y 1 — e+O.U!O931i ! +0.»00S16Oi 4 (1) 



\e(x)\ < 2.27 ■ 1CT 5 Vx > |e r (af)| < 1.21 • 1CT 4 Vx > 

Proof (only for absolute error). For < x < 4, see the figure above 0. 

For x > 4 let's consider that it is erf (4) = 0.99999998458... and erf(x) -> 1 

and erf is increasing, then 

(Vx>4) |1 -erf(x)\ < 10" 7 . (2) 

■ / -1.2735457x 2 -0.1487936x 4 , 

1 The figure represents the graph of 2.27- 10 -5 - \erf(x) - y 1 - e i+o.i48093ix 2 +o.ooo5i60*4 I 

made by the software Mathematica, showing that the quantity is positive. 

To verify that the graph do not intersect the x axis, you may make zooms in the subdomains 

[0.2,0.3], [0.75,0.85], [1.35,1.45] and [2.1,2.2], 



Let's T](x) our approximation of erf (2;) as in erf(x) = r](x) := y/T— e E ^ 

hmno- TTM -1.2735457a: 2 -0. 1487936a: 4 T+ • . 

ueing J^\-h) .— 1+0 . 148 093i.,;2 + 0.0005160x 4 • ±L IB ' 



(Va; > 4) < 3 = (^15 - l) 2 - 1 - 12 < • 4 2 - 0.795) 2 - 0.795 2 - 12 < 

= ~ ^°' 53 ) 2 ~ (^°' 53 ) 2 ~ 12 = 9~ x4 ~ °' 53x2 ~ 12 < 

< 0.14176a; 4 - 0.53x 2 - 12 = (0.148a; 4 - 0.00624)x 4 + (1.27a; 2 - 1.8)a; 2 - 12 



1.27a; 2 + 0.148a; 4 > 12 + 1.8a; 2 + 0.00624a; 4 



1.27a; 2 + 0.148a; 4 
12 



> 1 + 0.15a; 2 + 0.00052a; 4 



E(x) = - 



-1.2735457a; 2 - 0.1487936a; 4 



< 



-1.27a; 2 - 0.148a; 4 



+ 0.1480931a; 2 + 0.0005160a; 4 1 + 0.15a; 2 + 0.00052a; 4 



< -12 



,*(«) < e -u < £ -i2 + (e -i2 _ e - 24) = 2e - 



12 -24 



1 - (1 - e 



-12\2 



(1 - e- 12 ) 2 < 1 - e R ^ => 1 - e- 12 < y/l - 



< 1 - T)(x) = 1 - y/l - e R{x) 



< e 



-12 



(Va;>4) \l-r)(x)\<e 



(3) 



By (2) and Q it is 



(Vi > 4) |er/(at) - r/(a;)| < |1 - er/(a:)| + |1 - 7?(af)| < 10~ 7 + e~ 12 < 10 



□ 
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TABLE ERF2.27E-5 

Simply explicitly invertible approximations to 4 decimals of the error 
function erf, of complementary error function erfc, of the normal 
cumulative distribution function $, and of Q function, for x > 0. 
With majorizations of absolute and relative errors, and definitions. 



erf(x) 



=e _ * dt = 1 



erfc(x) <&(x) := 



1 



-oo V^n 



t = 1 - Q(x) 



A. 



erf(x) 



-1.2735457x2 -0.1 1S793G > 1 
\ _ g 1 + 0. 1480931x^+0. 0005160x 4 



|e(a;)| < 2.27- 10" 5 Vx > Q 
|e r (a;)| < 1.21 ■ 10~ 4 Vx > 



B. 



erfc(x) = 



-1.2735457x2 -0.1487936x 4 
\ — g 1 + 0. 1480931x^+0. 0005160x 4 



\e(x)\ < 2.27 ■ 10^ 5 Vx > 
s r (x)\ < 1% Vx 6 [0,6], 6 > 2.1588 



a $(x) ^ 



I 1 / -1.2735457^ 2 -0.074396Sic 4 
| y ]_ — g 2 + 0. 1480931^ 2 +0. 0002580x 4 



|e(x)j < 1.14- 10" 5 Vx > 
|e r (x)| < 1.78- 10~ 5 Vx > 



D. 



Q(x) = 



^ ^ / -1.2 7:15 lr>7 > - -0.0743968X 4 
V 1 — g2 + 0.1480931x 2 +0.0002580x 4 

2 2 



|e(x)| < 1.14 • 10" 5 Vx > 
|e r (x)| < 1% Vx e [0,6], b > 3.053 
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2 For x > 4.125 the approximation erf(x) = 1 has less absolute error. 
3 For x > 4.125 the approximation erfc(x) = has less absolute error. 
4 For x > 5.834 the approximation &(x) = 1 has less absolute error. 
5 For x > 5.834 the approximation Q(x) = has less absolute error. 
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